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We present the first measurements of identified iiadron production, azimutlial anisotropy, and 
pion interferometry from Au+Au collisions below the nominal injection energy at the Relativistic 
Heavy-Ion Collider (RHIC) facility. The data were collected using the large acceptance STAR 
detector at ^/snn = 9.2 GeV from a test run of the collider in the year 2008. Midrapidity results on 
multiplicity density (dN/dy) in rapidity (y), average transverse momentum ((pt)), particle ratios, 
elliptic flow, and HBT radii are consistent with the corresponding results at similar ^sjvat from fixed 
target experiments. Directed flow measurements are presented for both midrapidity and forward 
rapidity regions. Furthermore the collision centrality dependence of identified particle dN/dy, (pr), 
and particle ratios are discussed. These results also demonstrate that the capabilities of the STAR 
detector, although optimized for ^/sNN = 200 GeV, are suitable for the proposed QCD critical point 
search and exploration of the QCD phase diagram at RHIC. 

PACS numbers: 25.75.-q; 25.75.Dw; 24.85. +p; 25.75.Ld; 25.75.Gz 



I. INTRODUCTION 

Exploring the Quantum Chromodynamics (QCD) 
phase diagram is one of the goals of high energy 
heavy-ion collision experiments [1]. The QCD phase 
diagram is usually plotted as temperature (T) versus 
baryon chemical potential (/xb)- Assuming a thcr- 
malizcd system is reached in heavy-ion collisions, 
both of these quantities can be varied by changing 
the collision energy [2] . The phase diagram shows a 
possible transition from a high energy density and 
high temperature phase dominated by partonic de- 
grees of freedom, to a phase where the relevant de- 
grees of freedom arc hadronic [3]. Several obser- 
vations at the top RHIC energy, such as the sup- 
pression of high transverse momentum (px) hadron 
production in Au-|-Au collisions relative to p+p col- 
lisions [4], large elliptic flow {V2) for hadrons with 
light, as well as heavier strange valence quarks, and 
differences between baryon and meson V2 at inter- 
mediate pt for Au+Au collisions, have been asso- 
ciated with the existence of a phase with partonic 
degrees of freedom in the initial stages of heavy-ion 
collisions [1, 4, 5]. Lowering the collision energy and 
studying the energy dependence of these observables 
will allow us to search as a function of center of mass 
energy (-/s]vjv) or (T, /xb) for the onset of the tran- 
sition to a phase with partonic degrees of freedom 
at the early stage of the collision. 

Lattice QCD calculations indicate that the sys- 
tem produced at ^b = evolves through a rapid 
crossover in the quark-hadron phase transition [6]. 
Calculations from lattice QCD [7] and those from 
several QCD-based models [8] suggest that for col- 
lisions corresponding to large /ib, the transition is 
first order. The point in the (T, /xb) plane where 
the first order phase transition ends, is the QCD 
critical point [9] . Theoretical predictions of the loca- 
tion of this point on the phase diagram are subject 
to various ambiguities [10]. An experimental pro- 
gram for locating the QCD critical point through 
its signatures [10, 11] (e.g., long range fluctuations 
in event-by-event observables) is one of the exciting 



possibilities at the RHIC facility. These motivations 
form the basis of the proposal [12] by the experi- 
ments at RHIC to carry out a detailed program of 
exploring the phase diagram by varying the collision 
energy in high energy heavy-ion collisions. 

As an initial step to test the capabilities of the 
collider and experiments, a short run was conducted 
in the year 2008 at RHIC. The Au ions were col- 
lided at y/SNN — 9.2 GeV, which is below the in- 
jection energy of ^/snn = 19-6 GeV. At and be- 
low nominal injection energy, RHIC runs as a col- 
liding storage ring, further details of which can be 
found in Ref. [13]. The data taking period lasted 
for less than five hours at the Solenoidal Tracker at 
RHIC (STAR) experiment. This paper presents re- 
sults based on the analysis of this small data set and 
demonstrates the success of the test run in achiev- 
ing its objectives. The measurements shown here 
are the first step towards a detailed exploration of 
the QCD phase diagram at RHIC. 

The paper is organized as follows: The next sec- 
tion bricfiy presents the detectors used and details of 
the data analysis. In section HI, we present the re- 
sults including spectra. dN/dy, (pt) and particle 
ratios as a function of collision centrality and 
We also discuss results on directed fiow (wi), ehip- 
tic fiow (^2), and pion interferometry in this section. 
In section IV, we discuss the freeze-out conditions. 
Finally, in section V we summarize the results and 
provide a brief outlook of the upcoming Beam En- 
ergy Scan program at RHIC. 



II. EXPERIMENT AND DATA ANALYSIS 

A. STAR detector 

The results presented here are based on data taken 
at STAR [14] in Au-|-Au coUisions at y/sWW = 9.2 
GeV. This data set is taken with a minimum bias 
trigger. The trigger detectors used are the Beam- 
Beam Counter (BBC) and Vertex Position Detec- 
tor (VPD) [15]. The BBCs are scintillator annuh 
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mounted around the beam pipe beyond the east and 
west pole-tips of the STAR magnet at about 375 cm 
from the center of the nominal interaction region 
(IR). The inner tiles of the BBCs, with a pseudo- 
rapidity [t]) range of 3.8 < \r]\ < 5.2 and full az- 
imuthal coverage = 27r, are used to reconstruct 
the first-order event plane for the directed flow anal- 
ysis. The VPDs arc based on the conventional tech- 
nology of plastic scintillator read-out by photomul- 
tiplier tubes. They consist of two identical detector 
assemblies very close to the beam pipe, one on each 
side at a distance of \Vz\ = 5.6 m from the center 
of the IR. The main detector used to obtain the re- 
sults on pt spectra, yields, particle ratios, azimuthal 
anisotropy parameters, and pion interferometry for 
charged hadrons is the Time Projection Chamber 
(TPC) [16]. The TPC is the primary tracking de- 
vice at STAR. It is 4.2 m long and 4 m in diameter. 
Its acceptance covers ±1.8 units of pseudorapidity 
(77) and the full azimuthal angle. The sensitive vol- 
ume of the TPC contains PIO gas (10% methane, 
90% argon) regulated at 2 mbar above atmospheric 
pressure. The TPC data are used to determine parti- 
cle trajectories, momenta, and particle-type through 
ionization energy loss (dE/dx). STAR's solenoidal 
magnet field used for this low energy Au-|-Au test 
run was 0.5 T. In addition we present directed flow 
measurements from forward rapidities. These results 
used the data taken by the Forward Time Projection 
Chambers (FTPCs) [17]. There are two FTPCs lo- 
cated around the beam axis on both sides of the col- 
lision point. The sensitive medium is a gas mixture 
of equal parts Ar and CO2 by weight. The FTPCs 
detect charged particles in the pseudorapidity region 
2.5 < I77I < 4.0, with full azimuthal coverage. The 
details of the design and other characteristics of the 
STAR detectors can be found in Ref. [14]. 



B. Event selection 

The primary vertex for each minimum bias event 
is determined by finding the best point of common 
origin of the tracks measured in the TPC. The dis- 
tribution of the primary vertex position along the 
longitudinal beam direction (Vz) is shown in Fig. 1. 
The distribution is a broad Gaussian varying be- 
tween —200 and 200 cm, with a root mean square 
deviation of 89 cm. Only those events which have a 
Vz within 75 cm of the nominal collision point (cen- 
ter of the detector) are selected for the analysis, cor- 
responding to 57% of the total events recorded. This 
value is chosen by the trade-off between uniform de- 
tector performance within |?7| < 1.0 and sufhcicnt 
statistical significance of the measured observables. 
In order to reject events which involve interactions 
with the beam pipe and beam-gas interactions, the 
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FIG. 1: Event-by-event distribution of the z-position of 
the primary vertex (T4) in Au-|-Au collisions at y^sjvjv 
= 9.2 GeV. The vertical solid lines show the condition 
of IVzj < 75 cm for selected events. 



Vy where Vx 



event vertex radius (defined as ^JVx 
and Vy arc the vertex positions along the x and y 
directions) is required to be less than 2 cm. The 
Vx vs. Vy distribution is shown in Fig. 2. The circle 
with dotted lines corresponds to the event vertex ra- 
dius of 2 cm. A total of about 3000 events pass the 
selection criteria described above. 
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Au-|-Au collisions at -^sjvjv =9.2 GeV. The circle with 
dotted lines corresponds to a radius (= y^V^ + Vy) of 2 
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C. Centrality selection 



Centrality classes in Au+Au collisions at y^s^vtv 
= 9.2 GeV are defined using the number of charged 
particle tracks reconstructed in the main TPC over 
the full azimuth, pseudorapidity \r]\ < 0.5 and < 
75 cm. 

Figure 3 shows the uncorrected multiplicity distri- 
bution for charged tracks from the real data {Nj^^^^, 
open circles) and for those obtained from simulation 
(dashed histogram). Simulated multiplicity density 
is calculated using the two-component model [18] 
with the number of participants (A'part) and num- 
ber of collisions (A^coii) extracted from the Glauber 
Monte Carlo simulation as 



ch 



part 



-xN, 



coll 



(1) 



Here Upp is the average multiplicity in minimum bias 
p+p collisions and x is the fraction of the hard com- 
ponent. The inelastic cross-section for p+p colli- 
sions used in the Glauber Model simulations is 31.5 
mb [19]. The event- by-event multiplicity fluctua- 
tion has been taken into account by convoluting the 
Negative Binomial Distributions (NBD) for a given 
TVpart and iVcoU. The NBD distribution in multiplic- 



ity n has two parameters, Up 
as, 



and k, and is defined 



Pi 



{upp, fc; n) 



NBD {"-pp, 



r(n + k) 



ijipp/k)"' 



r(n + l)T{k) [upp/k + 1)«+'^' ' 

(2) 

where F is the Gamma function. The values k — 2.1 
and Upp = 1.12 are obtained by fitting the measured 
multiplicities with those from the simulation. The 
simulated multiplicity distribution is not sensitive to 
the k parameter. The distributions are found to be 
similar for varying k values such as k — 1.0, 1.6, and 
3.0. The fitting is performed for Nch > 17 in order 
to avoid the trigger inefficiency in peripheral colli- 
sions. The X value is fixed at 0.11 ± 0.03, obtained 
by extrapolating data from the PHOBOS collabo- 
ration [20]. The centrality is defined by calculating 
the fraction of the total cross-section obtained from 
the simulated multiplicity. 

Table I lists the centrality selection criteria for 
Au-I-Au collisions at ^snn = 9.2 GeV. We have 
divided the events into three centrality classes, 0- 
10%, 10-30%, and 30-60% of the total cross-section. 
The mean values of iVpart and A'coii have been eval- 
uated for these centrality bins and are given in 
Table I. Systematic uncertainties on (A'part) and 
(iVcoii) have been estimated by varying Upp and x 
in the two-component model as well as varying the 
input parameters in the Glauber Monte Carlo sim- 
ulation. The final errors on (A^pait) and {Nco\i) are 
the quadrature sum of these individual systematic 
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FIG. 3: Uncorrected charged particle multiplicity dis- 
tribution (open circles) measured in the TPC within 
|r;| < 0.5 in Au -I- An collisions at ^/sNN = 9.2 GeV. 
The dashed histogram represents the simulated multi- 
plicity distribution. The vertical dashed lines reflect the 
centrality selection criteria used in the paper. Errors are 
statistical only. 



TABLE I: Centrality selection, average number of partic- 
ipating nucleons ((A'^part)), and average number of binary 
collisions ((A'^coii)). 



% cross-section 




(A^part) 


(Acoll) 


0-10 


> 162 


317 ± 4 


716 ± 83 


10-30 


74-162 


202 ± 11 


395 ± 34 


30-60 


17-74 


88 ± 10 


133 ± 20 



errors. The results presented in this paper cover the 
collision centrality range of 0-60%. The results from 
more peripheral collisions are not presented due to 
large trigger inefficiencies in this test run. 



D. Track selection and particle identification 

Track selection criteria for the various analyses are 
presented in Table II. In order to avoid admixture 
of tracks from secondary vertices, a requirement is 
placed on the distance of closest approach (DCA) 
between each track and the event vertex. In order to 
prevent multiple counting of split tracks, a condition 
is placed on the number of track points (A^fit) used 
in the reconstruction of the track. Tracks can have 
a maximum of 45 hits in the TPC. 

To extract the pion yield in a given pT bin, we per- 
form an eight-Gaussian fit to the normalized dE/dx 
distributions of positively charged and negatively 
charged hadrons, simultaneously. The normalized 
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TABLE II: Track selection criteria for various analyses 
presented in this paper. 



Analysis 


DCA (cm) 






or y 


PT (GeV/c) 


PT spectra 


< 3 


> 20 


\y\ 


< 0.5 


> 0.1 


t'i(TPC) 


< 1 


> 20 


\v\ 


< 1.3 


0.15 - 2.0 


^;i(FTPC) 


< 1 


> 5 


2.5 < 


|r;| < 4.0 


0.15 - 2.0 


V2 


< 3 


> 15 


hi 


< 1.0 


0.1 " 2.0 


HBT 


< 3 


> 15 


\y\ 


< 0.5 
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FIG. 4: (Color online) dE/dx distribution for positively 
charged hadrons in the TPC, normalized by the expected 
pion dE/dx at 0.3 < pr < 0.4 GeV/c and \y\ < 0.5 
in Au+Au collisions at ^/snn ~ 9.2 GeV. The curves 
are Gaussian fits representing contributions from pions 
(dot-dashed, red), electrons (dashed, green), kaons (dot- 
dashed, blue), and protons (dotted, magenta). See text 
for details. Errors are statistical only. 



A similar procedure is followed to obtain yields for 
other pt ranges and for yields of kaons and protons. 
Further details of extracting raw yields of identified 
hadrons from normalized dE/dx distributions can 
be found in Ref. [23]. 

For the elliptic flow analysis of identified hadrons, 
the criteria of Incrj^l < 2 and \nap\ < 2 are used 
for extracting pion and proton V2. Since the mea- 
surements are carried out at low px (< 1.0 GeV/c), 
such an identification criterion is reasonable. For the 
pion interferometry analysis, the particle identifica- 
tion conditions are Itict^I < 2, \nap\ > 2, and \naK\ 
> 2, and the average transverse momentum (kr = 
{\piT + P2t|)/2) is required to fall in the range 150- 
250 MeV/c. 



E. Event plane for azimuthal anisotropy 

Azimuthal anisotropy can be quantified by study- 
ing the Fourier expansion of the azimuthal angle (</)) 
distribution of produced particles with respect to the 
reaction plane angle (^'i?) [24]. The various (order 
n) coefficients in this expansion are defined as: 



(cos 71 



(4) 



The angular brackets in the definition denote an av- 
erage over many particles and events. Directed fiow 
can be quantified by the first coefficient (vi) and el- 
liptic fiow by the second coefficient {V2), obtained 
using the above equation. 

In the azimuthal anisotropy analysis, vi and V2 are 
obtained from the following procedure. The event 
fiow vector (Qn) and the event plane angle (^'n) are 
defined by [24] 



Qn cos(n*,; 



E 



Wicos(n0i), (5) 



dE/dx in general is defined as 

\og{{dE/dx)/Bx) 
nax = , (3) 

where X is the particle type (e^,7r*,/C*, p, or p), 
Bx is the expected mean dE/dx of particle X, and 
ax is the dE/dx resolution of the TPC, which is 
a function of the track length in the TPC. The ex- 
pected mean dE/dx of particle X is calculated us- 
ing a Bichsel function for the energy loss in thin 
layers of PIO in the STAR TPC [16, 21]. Good 
agreement between the measurement and the cal- 
culation was demonstrated previously [22] . Figure 4 
shows a typical dE/dx distribution normalized to 
the pion dE/dx (referred to as the na^n- distribution) 
for charged hadrons with 0.3 < pr < 0.4 GeV/c and 
\y\ < 0.5. The counts under the Gaussian about naT^ 
~ give the yield of pions for a particular px range. 



ny 



tan 



E 



Wi sin(n</)^), (6) 



/n, 



(7) 



where sums go over all particles i used in the event 
plane calculation, and </>,; and Wi are the laboratory 
azimuthal angle and the weight for the i-th particle, 
respectively. Tracks used for the calculation of 
are excluded from the calculation of the event plane 
to remove self-correlation effects. 

Since finite multiplicity limits the angular resolu- 
tion of the reaction plane reconstruction, the 
has to be corrected for the event plane resolution by 



, ,obs 



(cos[n(*„ - *i^)]) 



(8) 
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to obtain the real w„ , where angular brackets denote 
an average over a large event sample. The event 
plane resolution is estimated from the correlation of 
the event planes of two sub-events. Assuming the 
pure flow correlations between the sub-events, the 
event plane resolution is given by 



{cos[ni^^ - ^iJn)]) = ^J{cos[n{^^~^B)]), (9) 

where A and B denote two subgroups of tracks. In 
this analysis, we use two random sub-events with 
equal numbers of particles. The full event plane res- 
olution is obtained from the resolution of the sub- 
events by 

(cos[n(*„ - = C{cos[n{^^ ~ *h)]>, (10) 

where C is a constant calculated from the known 
dependence of the resolution on multiplicity [24]. 

For the elliptic flow measurements presented in 
this paper, the TPC tracks are used to reconstruct 
the reaction plane [24] . The weights are taken to be 
the value of pr in GeV/c up to 2 GeV/c and then 
constant at 2.0 for pt > 2 GeV/c. Such weight val- 
ues are chosen as V2 increases with up to 2 GeV/c 
and then tends to saturates beyond p^ = 2 GeV/c. 
The variation of event plane resolution with collision 
centrality is shown in Fig. 5. The values of the reso- 
lution depend on the multiplicity and flow observed 
in the events. The resolution values are lower for 
y^J/vw = 9.2 GeV, compared to collisions at ^/snj\j 
= 200 GeV for similarly defined collision centrality 
classes [25]. A similar procedure for correcting the 
observed flow values with the resolution factor is fol- 
lowed for vi measurements. The vi results presented 
here are obtained using two different methods: the 
mixed harmonics and the standard methods. 

In the mixed harmonics method, vi is calcu- 
lated using mixed harmonics involving the second- 
harmonic event plane [25]. This method utilizes 
the large elliptic flow signal, and at the same time 
suppresses the non-flow contributions arising from 
the correlation of particles from the same harmon- 
ics. The method uses the second order event plane 
from the TPC (^'J^'") and the first order event plane 
from random sub-events in the FTPCs (vI'^'^'pc:! ^^j-^^^ 
^FTPC2)^ The average resolution for the event plane 
(as defined in Eqs. (9) and (10)) reconstructed from 
the TPC is 0.46 ± 0.03, while that reconstructed 
from the FTPCs is 0.41 ± 0.03, for 0-60% collision 
centrality. The mixed harmonics method is denoted 
by i;i{EPi,EP2} [25], as given below: 



7;i{EPi,EP2} 



(cos(( 

;cos(*^'^^^i 



+ ^FTPC _ 2^TPC)^ 

:*™^^^2^r^'))ite(*J^ ' 

(11) 



where the emission angle of the particle (0) is corre- 
lated with the vj/FTPC q£ ^j^^ random sub-event com- 
posed of tracks from both the FTPCs excluding that 
particle. 

In the standard method, the first-order event 
plane is reconstructed separately from the FTPC 
tracks (vi{EPi, FTPC}) and from the BBC hits 
(wi{EPi, BBC}). The event plane reconstructed 
from the detector on one side of the collision (east or 
west) is called a sub-event plane. A combination of 
the east and west sub-event plane vectors provides 
the full event plane. In the ■yi{EPi, FTPC} method, 
we used the event plane obtained from the full FTPC 
region to obtain the directed flow values measured 
in the TPC range {\t]\ < 1.3). A self-correlation 
arises if vi is obtained using particles from the same 
pseudorapidity region as is used for the event plane 
reconstruction. To avoid this self-correlation in the 
ui{EPi,FTPC} method, vi is obtained in the east 
FTPC (-4.2 < 7/ < -2.5) by using the event plane 
reconstructed in the west FTPC (2.5 < 77 < 4.2), 
and vice versa. In the '(;i{EPi,BBC} method, the 
event plane is obtained from the full BBC region (3.8 
< \r]\ < 5.2). Particles used for the estimation of vi 
with respect to the BBC full event plane cover the ij 
range up to 3.8, in order to avoid the self-correlation. 
The average resolution of the first order event plane 
for vi{EPi,FTPC} is 0.41 ± 0.03 for 0-60% cen- 
tral collisions while that for wijEPi, BBC} is 0.24 ± 
0.07. 




40-50 20-30 0-10 
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FIG. 5: Second order event plane resolution measured in 
the TPC as a function of collision centrality for Au+Au 
collisions at y^sjviv ~ 9.2 GeV. Errors are statistical only. 
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F. Correlation function in tt interferometry 

Experimentally, the two-particle correlation func- 
tion is obtained from the ratio, 



C(g,fc) = 



A{q, k) 
B{q, k) 



(12) 



We assume a spherical Gaussian source of 5 
fm [31, 32] for Au-|-Au collisions at ^/snn = 9.2 
GeV. The first term (1 — A) in Eq.(13) accounts 
for those pairs which do not interact or interfere. 
The second term represents those pairs where both 
Bose-Einstein effects and Coulomb interactions are 
present [32]. 



where A{q, k) is the distribution of particle pairs 
with relative momentum q = pi ^ P2 and average 
momentum k = {pi+p2)/2 from the same event, and 
B(q, k) is the corresponding distribution for pairs of 
particles taken from different events [26, 27]. The 
correlation function is normalized to unity at large 
q. In the mixed events, each particle in a given event 
is mixed with all particles (tt" for the results pre- 
sented in this paper) from other events, within a 
collection of 50 similar events. Similar events are se- 
lected within the centrality bin and further binned to 
have primary vertex z positions within 10 cm. With 
the availability of high statistics data and develop- 
ment of new techniques, it has become possible to 
have a three-dimensional decomposition of q [28-30] , 
providing better insight into the collision geometry. 

The relative momentum q can be decomposed ac- 
cording to the Bertsch-Pratt (also known as "out- 
side-long") convention [31]. The relative momentum 
q is decomposed into the variables along the beam 
direction (giong), parallel (^out) to the transverse mo- 
mentum of the pair kr = {pn +P2t)/2, and perpen- 
dicular ((?sidc) to (Jong and Qout- In addition to the 
correlation arising from quantum statistics of two 
identical particles, correlations can also arise from 
two-particle final state interactions. For identical 
pions, the effects of strong interactions are negligi- 
ble, but the long range Coulomb repulsion causes a 
suppression of the measured correlation function at 
small q. 

In this analysis, we follow the same procedure as 
was used in the previous analysis of Au-fAu colli- 
sions at -y/sATJv = 200 GeV [32] . For an azimuthally- 
integrated analysis at midrapidity in the longitudi- 
nal co-moving system (LCMS), the correlation func- 
tion in Eq. (12) can be decomposed as [33]: 

(^(^out, ^sido, fflong) — (1^^^) + 



I ' I ' I 



' I ' I ' I ' _ 

p.^ : 0.4 - 0.7 GeV/c 

• p Embedding _ 
o p Data 
□ p Data 




1 1.5 2 
DCA (cm) 



CO 



0.08 
0.06 
0.04 
0.02 




.'I'M 


1 ' 1 ' 1 ' 1 ' 1 . 


- 71 P-r 


0.4 - 0.7 GeV/c 


- • Embedding 


L o Data 


• • 






[ (b) 


^ cf - 




o - 




1,1,1,1,1" 



20 25 30 35 40 45 50 
Number of fit points 



FIG. 6: (a) Distribution of distance of closest approach 
of proton tracks to the primary vertex. The embedded 
tracks are compared to the ones in real data at 0.4 < 
Pt < 0.7 GeV/c at midrapidity in Au-|-Au collisions at 
^/sNN = 9.2 GeV. The DCA distribution of anti-protons 
in a similar kinematic range is also shown for comparison, 
(b) Comparison between the distributions of number of 
fit points for pions from embedding and from real data 
for 0.4 < Pt < 0.7 GeV/c at midrapidity in Au+Au 
collisions at ^/sNN = 9.2 GeV. 



A/CcOul(9i„v)(l + e-9o„t«Lt-9Lo«Le-?fo„,flU), 

(13) 

where i^coui is to a good approximation the squared 
nonsymmetrized Coulomb wave function integrated 
over a Gaussian source (corresponding to the LCMS 
Gaussian radii i?out, -Rsidc, ^long) [34]. Assuming 
particle identification is perfect and the source is 
purely chaotic, A represents the fraction of corre- 
lated pairs emitted from the collision. 



G. Correction factors for px spectra 

Two major correction factors for spectra ac- 
count for the detector acceptance and for the effi- 
ciency of reconstructing particle tracks. These are 
determined together by embedding the tracks sim- 
ulated using the GEANT [35] model of the STAR 
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FIG. 7: (a) Efficiency x acceptance for reconstructed 
pions, kaons, and protons in tlie TPC as a function of 
Pt at midrapidity in Au+Au collisions at ^/snn = 9.2 
GeV. (b) Percentage of pion background contribution 
estimated from HIJING+GEANT as a function of pr at 
midrapidity in Au+Au collisions at ^ysNN ~ 9.2 GeV. 
The contributions from different sources and the total 
background are shown separately. 



detector into real events at the raw data level. One 
important requirement is to have a match in the 
distributions of reconstructed embedded tracks and 
real data tracks for quantities reflecting track quality 
and used for track selection. Figures 6(a) and 6(b) 
show the comparisons of DCA (for protons) and A'fit 
(for pions) distributions, respectively, in the range 
0.4 < Pt < 0.7 GeV/c. Similar agreement as in 
Fig. 6 is observed between embedded tracks and real 
data in other measured px ranges for all the iden- 
tified hadrons presented in this paper. The ratio of 
the distribution of reconstructed and original Monte 
Carlo tracks as a function of pr gives the acceptance 
X efficiency correction factor as a function of px for 
the rapidity interval studied. The typical efficiency 
X acceptance factors in 0-60% central collisions for 
pions, kaons and protons at midrapidity {\y\ < 0.5) 



are shown in Fig. 7(a). The raw yields are corrected 
by these factors to obtain the final pT spectra. 

The STAR experiment has previously observed 
that proton yields had significant contamination 
from secondary protons, due to interactions of en- 
ergetic particles produced in collisions with detector 
materials. As these secondary protons are produced 
away from the primary interaction point, they ap- 
pear as a long tail in the DCA distribution of pro- 
tons. A comparison between shapes of DCA distri- 
butions of protons and anti-protons (which do not 
have such sources of background) was used in STAR 
to estimate the background contribution to the pro- 
ton yield [23, 36]. This feature was found to be 
more pronounced at lower pT- In this test run, it 
is observed that the DCA distribution for protons 
does not exhibit a long tail, and that for all the 
Pt ranges studied, its shape is similar to that for 
anti-protons (Fig. 6(a), distributions normalized to 
the same number of total counts). This lack of sec- 
ondary protons for Au+Au collisions at y^sJ^N = 9.2 
GeV could be due to the experimental configuration 
in the year 2008 with reduced amount of material in 
front of the STAR TPC, and due to the relatively 
small number of energetic particles produced in the 
interactions compared to collisions at higher ener- 
gies of ^/sNN = 62.4 and 200 GeV. No corrections 
for secondary proton background are applied for the 
present analysis at ^/snn ~ 9-2 GeV. 

The charged pion spectra are corrected for feed- 
down from weak-decays, muon contamination, and 
background pions produced in the detector mate- 
rials. These corrections are obtained from Monte 
Carlo simulations of HIJING events at ^^snn = 9.2 
GeV, with the STAR geometry for year 2008 and 
a realistic description of the detector response used 
in GEANT. The simulated events are reconstructed 
in the same way as the real data. The weak-decay 
daughter pions are mainly from Kg, and are iden- 
tified by the parent particle information accessible 
from the simulation. The muons from pion decay 
can be misidentified as primordial pions due to their 
similar masses. This contamination is obtained from 
Monte Carlo simulations by identifying the decay, 
which is accessible in the simulation. The weak- 
decay pion background and muon contamination ob- 
tained from the simulation are shown in Fig. 7(b), 
as a function of simulated pion pT for 0-60% central 
Au+Au colhsions at ^/snn = 9.2 GeV. The final 
pion spectra are corrected for this background ef- 
fect. 

The low momentum particles lose energy while 
traversing the detector material. The track recon- 
struction algorithm takes into account the Coulomb 
scattering and energy loss, assuming the pion mass 
for each particle. Therefore, a correction for the 
energy loss by heavier particles {K^, p and p) is 
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TABLE III: Sources of systematic errors on yields of var- 
ious produced hadrons. See section II H for more details. 



Hadron 


14 


cuts 


y 


correction 


PID 


extrapolation 


n 


3% 


3.2% 


2% 


5% 


5% 


3% 


K 


3% 


6.2% 


2% 


5% 


10% 


8% 


P 


3% 


5.4% 


10% 


5% 


4% 


15% 



needed. This correction is obtained from embedding 
Monte Carlo simulations. The largest change in re- 
constructed pt is found to be ^20 MeV/c at pT = 
200 MeV/c. For all results presented in this paper, 
the track pr is corrected for this energy loss effect. 



H. Systematic errors 

Systematic uncertainties on the spectra are esti- 
mated by varying cuts, and by assessing the purity of 
identified hadron sample from dE/dx measurements. 
In addition, the Gaussian fit ranges are varied to es- 
timate the systematic uncertainty on the extracted 
raw spectra. The point-to-point systematic errors 
are quoted in figure captions. The statistical and 
systematic errors are added in quadrature and plot- 
ted for most of the results unless otherwise specified. 
For integrated particle yields, extrapolating yields to 
unmeasured regions in pT is an additional source of 
systematic error. These are estimated by comparing 
the extrapolations using different fit functions to the 
Pt spectra. The detailed procedure is described in 
Ref. [23]. A summary of various sources of system- 
atic errors on the identified hadron yields for 0-60% 
centrality in Au-|-Au collisions at -sJsnn =9.2 GeV 
is given in Table III. The column titled in 
Table III represents the systematic errors obtained 
by varying the Vz range in the analysis, "cuts" lists 
systematic errors due to variation of DCA and iVfit 
cut values, "y" represents the systematic effect on 
yields due to a variation in rapidity range from ± 
0.5 to ± 0.2, "correction" includes the contribution 
to systematic errors from track reconstruction effi- 
ciency and acceptance estimates, "PID" represents 
the systematic errors associated with particle iden- 
tification (obtained by varying the dE/dx cuts and 
the range of Gaussian fits to normalized dE/dx dis- 
tributions), and "extrapolation" refers to the con- 
tribution of systematic errors from the different fit 
functions used for obtaining yields in unmeasured 
Pt ranges. In addition, the systematic error arising 
due to the pion background estimation (discussed in 
the previous subsection) is also calculated. It is of 
the order of 6%. The total systematic errors are of 
the order of 11%, 16%, and 20% for pion, kaon, and 
proton yields respectively. 



The systematic errors in the directed flow anal- 
ysis are obtained, (a) by exploiting the symmetry 
in the measurements for forward and backward re- 
gions with respect to 77 = 0, and (b) by compar- 
ing vi calculated from different methods with var- 
ious sensitivities to non-flow effects [25]. In (a), 
we average vi from the mixed harmonics method 
(ui{EPi, EP2}) and from the two standard methods 
(wi{EPi,FTPC} and ui{EPi, BBC}), as discussed 
in section II E, and take the difference between the 
magnitude of vi in the forward and backward region 
as the systematic error due to the unbalanced detec- 
tor response. We report an absolute error of ~7.8% 
in the FTPC range (2.5 < \r]\ < 4.2), and negligible 
error in the TPC range. In (b), we average the mag- 
nitude of vi in the forward and backward region, and 
take the maximum difference between results from 
the three methods as the systematic uncertainty. An 
absolute error of ^10% is found for the FTPC range, 
and ~50% relative error for the TPC range. The 
wijEPi, BBC} method in the TPC range (|r?| < 1.3) 
is more reliable compared to the other two methods. 
This is due to the large rj gap between the BBC and 
TPC detectors, which helps subtract the non-flow 
effect. The 77 gap between the BBC and the TPC is 
up to 2.6 units, while it is only 1.3 units between the 
FTPC and the TPC. The total absolute systematic 
error on the vi estimate is calculated as the quadra- 
ture sum of components (a) and (b), which is ~10% 
(absolute error), for the FTPC range and ^50% (rel- 
ative value) for the TPC range. 

The systematic errors on the elliptic flow param- 
eter are evaluated by varying the event vertex selec- 
tion along the beam direction, varying the DCA cut 
value, and by using the ij sub-event method. The 
total systematic error on V2 is approximately 10%. 

For the pion interferometry analysis, we study the 
following sources of systematic error: track merging, 
track splitting, size of the source used for Coulomb 
correction, particle identification, and pair accep- 
tance for pious of opposite charges. The estimated 
systematic errors are less than 10% for all radii in 
the 0-60% centrality bin for 150 < fcr < 250 MeV/c, 
similar to those in Refs. [31, 32]. 



III. RESULTS 
A. Transverse momentum spectra 

Figures 8(a) and 8(b) show the transverse momen- 
tum spectra for tt^ and p (p), respectively. Fig- 
ure 9(a) and 9(b) show the spectra for and K~ , 
respectively, in Au-f-Au collisions at ^ysNN — 9.2 
GeV. The results are shown for the collision central- 
ity classes of 0-10%, 10-30%, 30-60%, and 0-60%. 
The p spectrum is shown only for 0-60% central- 
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FIG. 8: Transverse momentum spectra for (a) charged 
pions and (b) protons at midrapidity < 0.5) in 
Au+Au collisions at y^ijvjv = 9.2 GeV for various cen- 
tralities. The distributions for anti-protons were mea- 
sured in this limited statistics data only for 0-60% cen- 
trality. The anti-proton yield shown in the figure is mul- 
tiplied by a factor of 10. The errors shown are statistical 
and systematic errors (discussed in section II H) added 
in quadrature. 



FIG. 9: Transverse momentum spectra for (a) positive 
kaons and (b) negative kaons at midrapidity < 0.5) 
in Au-f Au collisions at ^/sNN = 9.2 GeV for various cen- 
tralities. The errors shown are statistical and systematic 
errors (discussed in section II H) added in quadrature. 



B. Centrality dependence of particle 
production 



ity and the yields are multiplied by a factor of 10 for 
visibility. The inverse slopes of the identified hadron 
spectra follow the order tt < K < p. An exponential 
fit to the pt spectra of 7r+, K^, and p yields inverse 
slopes of 180 ± 7 McV, 360 ± 7 McV and 616 ± 11 
McV respectively. The errors on the inverse slopes 
are statistical. The spectra can be further character- 
ized by looking at the dN/dy and (pt) or (mr) — m 
for the produced hadrons, where m is the mass of 
the hadron and ttit ~ \/ + p^ is its transverse 
mass. Those observables are discussed in the follow- 
ing sections. 



Figures 10 and 11 show the comparison of colli- 
sion centrality dependence of dN/dy of 7r+, K^, and 
p, normalized by (A'part)/2, between new results at 
^ysNN = 9.2 GeV and previously published results 
at y/sWN = 62.4 and 200 GeV from the STAR ex- 
periment [4, 23, 37]. The yields of charged pions 
and kaons decrease with decreasing collision energy. 
The collision centrality dependence within the lim- 
ited centrality region studied for the new results is 
similar to that at higher beam energies. For pro- 
tons, the yield is larger in central Au-|-Au collisions 
at y/sNN = 9.2 GeV compared to corresponding re- 
sults at ^i]^ = 62.4 and 200 GcV [4, 23, 37]. For 
the most peripheral collisions, the yields are compa- 
rable within errors to corresponding yields at higher 
beam energies. The increase in proton yield with 
the increasing collision centrality is due to large net- 
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FIG. 10: dN/dy of (a) tt^ and (b) p, normalized by 
{A''part)/2, for Au+Au collisions at ^ysNN = 9.2 GeV, 
plotted as a function of (A^part). The lower energy re- 
sults are compared to corresponding results for Au+Au 
collisions at y/s]^ = 62.4 and 200 GeV [23, 37]. The 
errors shown are the quadrature sum of statistical and 
systematic uncertainties. The systematic errors on pion 
and proton yields for ^/sNN = 9.2 GeV data are ~ 12% 
and ~ 20%, respectively, for all the collision centralities 
studied. 



FIG. 11: dN/dy of (a) K+ and (b) K~ , normalized by 
{A''part)/2 for Au+Au collisions at ^/snn = 9.2 GeV, 
plotted as a function of (Apart). The lower energy re- 
sults are compared to corresponding results for Au+Au 
collisions at yi]^ = 62.4 and 200 GeV [4, 23, 37]. The 
errors shown are the quadrature sum of statistical and 
systematic uncertainties. The systematic errors on Jf"*" 
and K~ yields for ^/sNN = 9.2 GeV data are similar, 
about 18% for all the collision centralities studied. 



TABLE IV: Centrality dependence of dNch/drj at midra- 
pidity in Au+Au collisions at ^^snn ~ 9.2 GeV. 



% cross section 


dNch/drj 


stat. error 


sys. error 


0-10 


229 


25 


62 


10-30 


133 


15 


36 


30-60 


48 


5 


13 



proton [p — p) density at midrapidity in the lower 
collision energies. 

The inclusive dNch/drj (sum of contributions from 
TT^, , and p (p) found by redoing the analysis 
binned in rj instead of rapidity) at midrapidity for 
various collision centralities arc given in Table IV 
along with the statistical and systematic errors for 
Au+Au collisions at y^SNN = 9.2 GeV. 



Figures 12 and 13 show the comparison of (px) as 
a function of (A'part) for 7r+, K~^, andp from Au+Au 
collisions at ^Sjvjv = 9.2 GeV with those from colli- 
sions at = 62.4 and 200 GeV [4, 23, 37]. For 
the collision centralities studied, the dependencies of 
(pt) on (A^part) at ^/snn = 9.2 GeV are similar to 
those at y^sjvlv = 62.4 and 200 GeV. An increase 
in (pt) with increasing hadron mass is observed at 
^/sNN =9.2 GcV. A similar dependence is also ob- 
served for ^/SNN = 62.4 and 200 GeV. However, 
the differences in (px) between protons and kaons 
are much smaller compared to the observations at 
higher beam energies. The mass dependence of (pt) 
reflects collective expansion in the radial direction. 
The smaller difference between (pt) of protons and 
kaons at ^/SNN =9.2 GeV indicates that the average 
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FIG. 12: (pt) for 7r+, , and p plotted as a function 
of (A'part) for Au+Au collisions at ^/sNN = 9.2 GeV 
and compared to corresponding results at y^sjviv = 62.4 
GeV [4, 23, 37]. The errors shown are the quadrature 
sum of statistical and systematic uncertainties. The sys- 
tematic errors for pions, kaons, and protons for y^ijviv 
= 9.2 GeV are ~ 12%, 18%, and 21% respectively, and 
similar for all the collision centralities studied. 
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FIG. 13: (pt) for tt"*", K'^ and p plotted as a function 
of (A'part) for Au+Au collisions at ^/snn = 9.2 GeV 
and compared to corresponding results at ^snn ~ 200 
GeV [4, 23, 37]. The errors shown are the quadrature 
sum of statistical and systematic uncertainties. The sys- 
tematic errors for pions, kaons, and protons for ^/snn 
= 9.2 GeV are ~ 12%, 18%, and 21% respectively, and 
similar for all the collision centralities studied. 



collective velocity in the radial direction is smaller 
at that energy. 

Figures 14 and 15 show the various particle ratios 
{K~/K^, K~/Tr~, p/n^, and K^/n^) as a func- 
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FIG. 14: Variation of (a) / K+ and (b) A'^/tt^ ratios 
as a function of (Apart) for Au+Au collisions at ^/snn = 
9.2 GeV. For comparison we also show the correspond- 
ing results from Au+Au collisions at ^/snn = 62.4 and 
200 GeV [23, 37]. The errors shown are the quadrature 
sum of statistical and systematic uncertainties. The sys- 
tematic errors for K~ /K^ and K~ /■k~ for ^snn = 9.2 
GeV data are ~ 25% and 22%, respectively, and similar 
for all the collision centralities studied. 



tion of collision centrality expressed as (Npai-t) for 
Au+Au collisions at ^/snn = 9.2 GeV. Correspond- 
ing results from Au+Au collisions at ^/snn = 62.4 
and 200 GeV [4, 23, 37] are also shown. The 7r-/7r+ 
ratio is close to unity and is not shown. Due to low 
event statistics and the low yield of anti-protons, the 
centrality dependence of the p/p ratio for y/s^N = 
9.2 GeV collisions could not be extracted. 

The / and A'^/tt^ ratios are lower at 
y/sJjN — 9.2 GeV compared to those at ^/Tnn = 
62.4 and 200 GeV. In the case oiK+/TT+, there is less 
variation between 9.2 GeV and the highest RHIC 
energies than in case of the other particle ratios dis- 
cussed above. This reflects an interplay between the 
decreasing importance of associated production and 
an increasing contribution from pair production of 
kaons with increasing collision energy. Associated 
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FIG. 15: Variation of (a) p/tt^ and (b) /-k^ ratios as 
a function of (A'^part) for Au+Au collisions at ^/snn = 
9.2 GeV. For comparison we also show the correspond- 
ing results from Au+Au collisions at ^/snn = 62.4 and 
200 GeV [4, 23, 37]. The errors shown are the quadra- 
ture sum of statistical and systematic uncertainties. The 
systematic errors for p/tv^ and K'^ /tt'^ for ^sjvjv ~ 9.2 
GeV data are ~ 25% and 22%, respectively, and similar 
for all the collision centralities studied. 



production refers to reactions such as NN — >■ KYN 
and TriV KY, where iV is a nucleon and Y a hy- 
peron. The p/n^ ratio is larger at ^/snn = 9.2 GeV 
than at ^^snn = 62.4 and 200 GeV for all collision 
centralities studied. As discussed above, this is a 
consequence of higher net-proton density at midra- 
pidity for the collisions at y^s^N = 9.2 GeV com- 
pared to those at ^/snn = 62.4 and 200 GeV. 



C. Energy dependence of particle production 

Figure 16 shows the dNch/df] at midrapidity nor- 
malized by (A'pait)/2 as a function of y^SNN- The 
result from y^jvw = 9.2 GeV is in agreement with 
the general energy dependence trend observed at the 
AGS [381, SPS [391, and RHIC [23, 401. The result 
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FIG. 16: The midrapidity dNch/drj normalized by 
(A'^part)/2 as a function of y^ijvjv. Au+Au collisions at 
y^SNN = 9.2 GeV are compared to previous results from 
AGS [38], SPS [39], and RHIC [23, 40]. The error shown 
are the quadrature sum of statistical and systematic un- 
certainties. The systematic error for ^/snn ~ 9.2 GeV 
data is 39%. 



at 9.2 GeV has a value close to that obtained at 
a similar energy {^/snn = 8.8 GeV) by the NA49 
experiment at SPS [39]. Figures 17(a) and 18(a) 
show dN/dy normalized by {Np^rtj/^ for tt^ and 
K^, respectively, in 0-10% central Au+Au collisions 
at y^SNN = 9.2 GeV, compared to previous results 
at AGS [38], SPS [39], and RHIC [23]. Within er- 
rors, the yields are consistent with previous results 
at similar ^/snn. Figures 17(b) and 18(b) show the 
(ttit) — m for tt^ and , respectively, in 0-10% 
central Au+Au colhsions at ^/sNN = 9.2 GeV. The 
results are also compared to previous measurements 
at various energies. The results from Au+Au colli- 
sions at y'sNN = 9.2 GeV are consistent with cor- 
responding measurements at SPS energies at similar 
y/sJjN. Both dN/dy and (my) — rn are obtained 
using data in the measured pt ranges and extrap- 
olations assuming certain functional forms for the 
unmeasured pT ranges, as discussed in section V.B 
of our previous paper [23] . For the present midra- 
pidity measurements, the percentage contribution to 
the yields from extrapolation are about 20% for tt^ , 
50% for K^, and 25% for p. 

The (my) — m values increase with ^/snn at 
lower AGS energies, stay independent of ^/snn at 
the SPS and RHIC 9.2 GeV collisions, then tend 



to rise further with increasing ^/sNN at the higher 
beam energies at RHIC. For a thermodynamic sys- 
tem, (mr) — rn can be an approximate representa- 
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FIG. 17: (a) dN/dy normalized by (iVpart)/2 and (b) 
(mr) — m of tt^, in 0-10% central Au+Au collisions for 
y/SNN ~ 9.2 GeV compared to previous results from 
AGS [38], SPS [39], and RHIC [23]. The errors shown 
are the quadrature sum of statistical and systematic un- 
certainties. 



FIG. 18: (a) dN/dy normalized by (iVpart)/2 (b) (mr) - 
m of K , in 0-10% central Au+Au collisions for 
^/sNN = 9.2 GeV compared to previous results from 
AGS [38], SPS [39], and RHIC [23]. The errors shown 
are the quadrature sum of statistical and systematic un- 
certainties. 



tion of the temperature of the system, and dN/dy 
(X ln(yj/v7v) may represent its entropy. In such a 
scenario, the observations could reflect the charac- 
teristic signature of a first order phase transition, 
as proposed by Van Hove [41]. Then the constant 
value of (tot) — to vs. ^/snn around 9.2 GcV has 
one possible interpretation in terms of formation of 
a mixed phase of a QGP and hadrons during the 
evolution of the heavy-ion system. However, there 
could be several other effects to which (tot) — to 
is sensitive, which also need to be understood for 
proper interpretation of the data [42]. The energy 
dependencies of the proton dN/dy and (tot) —to are 
not discussed in this paper, as the STAR results are 
presented without correction for feed down contri- 
butions. The low event statistics in the present data 
does not allow us to obtain feed-down corrections 
from the data itself. All results presented in this pa- 
per arc from inclusive protons and anti-protons as in 
our previous paper at higher energies at RHIC [23]. 



Figures 19(a) and 19(b) show the collision energy 
dependence of the particle ratios 7r^/7r+ andp/p, re- 
spectively, in central heavy-ion collisions. Similarly, 
Figs. 20(a) and 20(b) show the ratios of A'^//v+ and 
K/n^ respectively. The new results from Au+Au 
collisions at ^Jsnn — 9.2 GeV follow the ^Jsnn 
trend established by previous measurements. The 
PT-hitegrated 7r~/7r+ ratio at ^Jsmn = 9.2 GeV is 
1.08 ± 0.04 (stat.) ± 0.16 (sys.). Those at lower 
beam energies have values much larger than unity, 
which could be due to significant contributions from 
resonance decays (such as from A baryons). The 
value of the p/p ratio at ^/sNN = 9.2 GeV is 0.010 ± 
0.001 (stat.) ± 0.003 (sys.) indicating large values of 
net-protons. The p/p ratio increases with increasing 
collision energy and approaches unity for top RHIC 
energies. This indicates that at higher beam energies 
the p (p) production at midrapidity is dominated by 
pair production. The K~/K^ ratio at ^/snn = 
9.2 GeV is 0.38 ± 0.05 (stat.) ± 0.09 (sys.), indi- 
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FIG. 19: (a) tv^ /tv^ and (b) p/p ratios at midrapid- 
ity (I y l< 0.5) for central 0-10% Au+Au collisions at 
y/SNN ~ 9.2 GeV compared to previous results from 
AGS [38], SPS [39], and RHIC [23]. The errors shown 
are the quadrature sum of statistical and systematic un- 
certainties. The systematic errors on -k~ /n^ and p/p for 
y'sNN = 9.2 GeV data are 15% and 27%, respectively. 



FIG. 20: (a) K~ / and (b) K/tt ratios at midrapid- 
ity (I y l< 0.5) for central 0-10% Au+Au collisions at 
y/SNN = 9.2 GeV compared to previous results from 
AGS [38], SPS [39], and RHIC [23]. The errors shown 
are the quadrature sum of statistical and systematic un- 
certainties. The systematic errors on K~ / K'^ and K/-k 
for Y^sjviv = 9.2 GeV data are 23% and 19%, respectively. 



eating a significant contribution to kaon production 
from associated production at lower collision ener- 
gies. With increasing ^Jsnn, the /K^ ratio ap- 
proaches unity, indicating dominance of kaon pair 
production. The K/tt ratio is of interest, as it ex- 
presses the enhancement of strangeness production 
relative to non-strange hadrons in heavy-ion colli- 
sions compared to p -\- p collisions. The increase in 
ii'+/7r+ ratio with beam energies up to i/snn = 
7.7 GeV at SPS and the subsequent decrease and 
possible saturation with increasing beam energies 
has been a subject of intense theoretical debate re- 
cently [39, 43]. The discussions mainly focus on the 
question of the relevant degrees of freedom that are 
necessary to explain the energy dependence of the 
K/tt ratio. Our new results from Au+Au collisions 
at y^SNN = 9.2 GeV with only about 3000 events 
(hence with large errors) are found to be consistent 
with the previously observed energy dependence. 



D. Azimuthal anisotropy 

The study of collective flow in relativistic nuclear 
collisions could provide insights into the equation 
of state (EOS) of the matter created by the col- 
lisions. As discussed earlier, there are two types 
of azimuthal anisotropy that are widely studied in 
heavy-ion collisions, directed flow (vi) and elliptic 
flow (^2). Directed flow measurements at forward 
rapidities describe the "side-splash" motion of the 
collision products. Hence, it is an important tool 
to probe the dynamics of the system at forward ra- 
pidities [44]. Since vi is generated very early in the 
evolution of heavy-ion collisions, it probes the on- 
set of bulk collective dynamics. The shape of vi vs. 
rapidity around midrapidity is suggested as a signa- 
ture of a first order phase transition [45]. On the 
other hand, the characterization of the elliptic flow 
of produced particles by their azimuthal anisotropy 
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has proven to be one of the most successful probes of 
the dynamics in Au+Au colhsions at RHIC [25, 46- 
52]. Elhptic flow provides the possibiHty to gain 
information about the degree of thermahzation of 
the hot, dense medium. Studying its dependence 
on system size, number of constituent quarks, trans- 
verse momentum, and transverse mass, is crucial to 
the understanding of the properties of the produced 
matter. 
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FIG. 21: (Color online) Charged hadron vi vs. 77 from 
the 0-60% collision centrality Au+Au collisions at y^sjviv 
— 9.2 GeV. The errors shown are statistical. Systematic 
errors are discussed in section II H. The solid star sym- 
bols are the results obtained from the mixed harmonic 
method, while the open star and open plus symbols rep- 
resent results from the standard methods (see text for 
details). The results are compared to vi from 30-60% 
collision centrality Au-f Au collisions at y^sjvjv ~ 62.4 
and 200 GeV [53]. For comparison, vi for charged pions 
for the 0-60% collision centrality from Pb+Pb collisions 
at y^SNN = 8.8 GeV are also shown [54]. 



Figure 21 shows charged hadron vi results in 
Au+Au collisions for the 0-60% collision central- 
ity at ^/snn = 9.2 GeV, compared to correspond- 
ing results from 30-60% central Au+Au collisions at 
y/sW]^ = 62 A and 200 GeV [53]. The pr range of 
this study is 0.15-2.0 GeV/c. The vi results from 
Au+Au collisions at ^/sNN = 9.2 GeV are shown for 
the three different methods, as described in section 
II E. The results from the three methods are con- 
sistent within the error bars. These results are also 
compared with vi for charged pions in Pb+Pb colli- 
sions at ^/SNN = 8.8 GeV measured by NA49 [54]. 
At midrapidity, all the results have comparable val- 
ues. At forward rapidity {\ri\ > 2), the trend of vi 
for higher ^/sNN (62.4 and 200 GeV) appears to be 
different from that for y^JWW = 9.2 GeV. This can 
be explained by contributions from spectator pro- 
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FIG. 22: (Color online) Same as Fig. 21, but plotted as 
a function of ry/j/bcam. 



tons to the directed flow signal at large [77[. The 
beam rapidities (yboam) for ^/s^n = 9.2, 62.4, and 
200 GeV are 2.3, 4.2, and 5.4 respectively. With -q 
divided by the respective ybeam values for the beam 
energies (Fig. 22), all the vi values follow a common 
trend for the measured [?7[/?/beam < 1 range. 

Figure 23 shows V2{pt) for charged hadrons, pi- 
ons, and protons in ^/snn = 9.2 GeV collisions. For 
comparison, we show pion V2 results from NA49 [54] 
at similar y^SNN- Within the statistical errors, there 
is good agreement between results from the two ex- 
periments. At top RHIC energies, V2 at low px 
shows a characteristic scaling with particle mass [5] 
that is consistent with hydrodynamic behavior; how- 
ever, the available statistics in the current analysis 
are insufficient to extend this study to 9.2 GeV. The 
small number of events also precludes the extension 
of the measurements to larger px values, to study the 
number of constituent quark scaling of V2 observed 
at yijviv = 200 GeV. Figure 24 shows the elliptic 
flow parameter at ^/snn = 9.2 GeV compared to 
other beam energies [46, 50, 54, 56-58]. The STAR 
data at ^/snn = 9.2 GeV, denoted by the star sym- 
bol, follow the observed y'ijviv dependence. 



E. Pion interferometry 

Information about the space-time structure of the 
emitting source can be extracted with intensity in- 
terferometry techniques [59]. The primary goal of 
pion interferometry, performed at midrapidity and 
at low transverse momentum, is to study the space- 
time size of the emitting source and freeze-out pro- 
cesses of the dynamically evolving collision fireball. 
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FIG. 23: V2 as a function of pr for charged hadrons (solid 
triangles), tt (solid circles), and p (solid squares) in 0- 
60% Au+Au collisions at ^/sNN = 9.2 GeV. The error 
bars include only statistical uncertainties for y^sjviv = 
9.2 GeV data. The corresponding systematic error is 
discussed in section II H. For comparison, V2{pt) results 
for TT (open circles) from NA49 [54] in 0-43.5% Pb + Pb 
collisions at ^/sNN = 8.8 GeV, are also shown. 
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Au+Au collisions. Only statistical errors are shown, 
results of STAR charged hadron V2 [55] are compared to 
those measured by E877 [56], NA49 [54], PHENIX [57], 
and PHOBOS [46, 50, 58]. 



IV. FREEZE-OUT PARAMETERS AND 
PHASE DIAGRAM 



TABLE V: The HBT parameters for 0-60% central 
events and 150 < kr < 250 MeV/c. 



A 


-Rout (fm) 


-Rside (fm) 


-Riong (fm) 


0.36 ± 0.08 


5.05 ± 0.96 


3.52 ± 0.56 


3.25 ± 0.86 



The 3-diniensional correlation functions are fitted 
with Eq. (13), where Ri is the homogeneity length in 
the i direction [28-30]. Projections of the fit to the 
correlation function of the 0-60% most central colli- 
sions, weighted according to the mixed-pair back- 
ground, are shown in Fig. 25. The three panels 
show the projections of the 3-dimensional correla- 
tion function onto the (?out, (Zsidc, and giong axes. 
The curves show Bowler- Sinyukov fits [34] to the 
Coulomb-corrected correlation function. Table V 
lists the HBT parameters obtained from the fits 
along with statistical errors. 

The radius parameter -Rgido has the most direct 
correlation with the source geometry, whereas -Rout 
encodes both geometry and time scale. Hydrody- 
namic calculations with a first order phase transition 
predict a ratio of -Rout /-Rsido larger than unity. Our 
measurements indicate the ratio -Rout/-Rsidc = 1.4 ± 
0.4. 



The measured hadron spectra reflect the proper- 
ties of the bulk matter at kinetic freeze-out, after 
elastic collisions among the hadrons have ceased. 
More direct information on the earlier stages can 
be deduced from the integrated yields of the differ- 
ent hadron species, which change only via inelastic 
collisions. The point in time at which these inelas- 
tic collisions cease is referred to as chemical freeze- 
out, which takes place before kinetic freeze-out. The 
transverse momentum distributions of the different 
particles contain two components, one random and 
one collective. The random component can be iden- 
tified as the one that depends on the temperature of 
the system at kinetic freeze-out (Tkin). The collec- 
tive component, which arises from the matter den- 
sity gradient from the center to the boundary of the 
fireball created in high energy nuclear collisions, is 
generated by collective flow in the transverse direc- 
tion, and is characterized by its velocity (3t- 

Assuming that the system attains thermal equilib- 
rium, the blast-wave (BW) formulation [60] can be 
used to extract Tkin and (/3t)- The transverse flow 
velocity of a particle at a distance r from the center 
of the emission source, as a function of the surface 
velocity (/3s) of the expanding cylinder, is parame- 
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terized as Prir) 
fitting the data, 
trum is then 

dN 



prdpi 



= /3s(r/i?)", where n is found by 
The transverse momentum spec- 



, r /pTsinhp(r) 
r ar niTio 





xKi 



niT coshp(r) 

Tkin 



(14) 



where Iq and Ki are modified Bessel functions and 
p(r) = tanh~^ /3t(?')- Simultaneous fits to tlie px 
distributions of tt, K, and p at midrapidity for cen- 
tral 0-10% Au-f Au colhsions at ^/snn = 9.2 GeV 
are shown in Fig. 26. The extracted parameters are 
Tkin = 105 ± 10 (stat.) ± 16 (sys.) MeV, (/3t) = 
0.46c ± 0.01c (stat.) ± 0.04c (sys.), and ti = 0.9 ± 
6.4 (stat.) ± 6.4 (sys.) with xVndf = 15/17. The 
parameter n is poorly constrained by the fits in this 
low event statistical data set. The parameter val- 
ues do not change within the quoted errors for other 
centrality ranges. Only statistical errors are used 
for obtaining the fit parameters shown in the figure. 
Inclusion of systematic errors gives similar values of 
Tkin and (/3t). Similar studies have been done for 
other higher energy collisions at RHIC [23]. 



FIG. 25: Projections of the three-dimensional correla- 
tion function and corresponding Bowler-Sinyukov [34] 
fits (lines) for negative pions from the 0-60% central 
Au+Au events and kr = [150, 250] MeV/c. 
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FIG. 26: Midrapidity transverse momentum distribu- 
tions of pions, kaons, and protons (no feed-down correc- 
tion) for 0-10% most central Au-|-Au collisions at ^/snn 
— 9.2 GeV, fitted to blast- wave model calculations [60]. 
The extracted kinetic freeze-out parameters are Tkin = 
105 ± 10 (stat.) ± 16 (sys.) MeV and (/3t) = 0.46c ± 
0.01c (stat.) ± 0.04c (sys.). 
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FIG. 27: Midrapidity particle ratios for 0-10% most cen- 
tral Au-fAu coUisions at ^/snn = 9.2 GeV, fitted to 
thermal model calculations. See text for details. The 
extracted chemical freeze-out temperature is Teh = 151 
± 2 (stat.) ± 7 (sys.) MeV and baryon chemical poten- 
tial is fiB = 354 ± 7 (stat.) ± 30 (sys.) MeV. 

Within a statistical model in thermodynamical 
equilibrium, the particle abundance in a system of 
volume V can be given by 



9i 



(2^) 



3 7|' 



exp 



(^ 



± 1 



-d^p, 



(15) 
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where Ni is the abundance of particle species i, gi 
is the spin degeneracy, Bi and Si are the baryon 
number and strangeness number, respectively, Ei is 
the particle energy, and the integral is taken over 
all momentum space [23]. The model parameters 
are the chemical freeze-out temperature (Teh), the 
baryon {(Ib) and strangeness (fJ-s) chemical poten- 
tials, and the ad hoc strangeness suppression factor 
(75). Measured particle ratios are used to constrain 
the values of temperature (Tch) and baryon chemical 
potential (hb) at chemical freeze-out using the sta- 
tistical model assumption that the system is in ther- 
mal and chemical equilibrium at that stage. Fits are 
performed to the various ratios for midrapidity cen- 
tral 0-10% Au4-Au collisions at ^/sWn = 9.2 GeV 
using such a model, and are shown in Fig. 27. The 
analysis is done within the framework of a statisti- 
cal model as discussed in Ref. [61]. This model has 
been used to extract chemical freeze-out parameters 
at higher RHIC energies [23] . The extracted param- 
eter values are Tch = 151 ± 2 (stat.) ± 7 (sys.) 
MeV, fiB = 354 ± 7 (stat.) ± 30 (sys.) MeV, us = 
25 ± 9 (stat.) ± 14 (sys.) MeV, and 7s = 0.9 ± 0.7 
(stat.) ± 0.1 (sys.) for 9.2 GeV data. These values 
are very close to those extracted from the measure- 
ments at SPS for similar ^/snn [62]. Only statistical 
errors on the particle production ratios are used for 
obtaining the fit parameters. Inclusion of systematic 
errors gives similar values of Tch and /ib- 

Figure 28 shows the temperatures at various 
stages in heavy-ion collisions as a function of fiB 
(at different ^/snn)- The /iB values shown are es- 
timated at chemical freeze-out. The initial temper- 
atures (Tinitiai) achieved at top RHIC and SPS en- 
ergies are obtained from models [63] that explain 
the direct photon measurements from the PHENIX 
experiment at RHIC [64] and from the WA98 ex- 
periment at SPS [65]. From these models, which 
assume that thermalization is achieved in the colli- 
sions within a time between 0.1-1.2 fm/c, the Tinitiai 
extracted is greater than 300 MeV at RHIC and 
greater than 200 MeV at SPS. The Tch and Tki„ val- 
ues extracted from particle ratios and pr spectra of 
various hadrons, respectively, from models assuming 
thermodynamical equilibrium are also shown. The 
values for ^/snn ~ 9.2 GeV are from the data pre- 
sented in this paper. The values at other ^/snn are 
from Ref. [23] and references therein. It is inter- 
esting to observe that Tch and Tkin values approach 
each other in the high /iB regime. A few recent pre- 
dictions from lattice QCD calculations [10] are also 
shown in Fig. 28. Several lattice QCD calculations 
indicate that the partonic to hadronic phase transi- 
tion occurs around Tc ~ 170-190 MeV [66]. These 
calculations also suggest that the phase transition at 
/iB = is a cross-over [6]. Most QCD-based model 
calculations [3, 7] suggest that the phase transition 
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FIG. 28: Temperature vs. baryon chemical potential 
(/ie) from heavy-ion collisions at various y^ijvjv [23]. 
The values shown are estimated at chemical freeze- 
out. The kinetic and chemical freeze-out parameters, 
extracted using models assuming thermal and chemi- 
cal equilibrium from midrapidity measurement in cen- 
tral 0-10% Au-|-Au collisions at ^/snn = 9.2 GeV, are 
shown as star symbols. The range of critical tempera- 
tures (Tc) [66] of the cross-over quark-hadron phase tran- 
sition at /iB = [6] and the QCD critical point from 
two different calculations [9] from lattice QCD are also 
indicated. Model-based estimates of the range of ini- 
tial temperature (Tinitiai) achieved in heavy-ion collisions 
based in part on direct photon data at top RHIC [64] and 
SPS [65] energies are also shown. The range of /ib to be 
scanned in the upcoming RHIC critical point search and 
Beam Energy Scan program corresponding to y^sjvjv = 
5.5 to 39 GeV is indicated by horizontal arrows near the 
/IB axis [12]. 



at large fj,B is of first-order. Two estimates of the 
QCD critical point [9] in the T — /iB plane taking 
Tc = 176 MeV are shown in Fig. 28. The region 
planned to be explored in the critical point search 
program at RHIC is shown in Fig. 28. 



V. SUMMARY AND OUTLOOK 

We have presented measurements of identified 
particle production, azimuthal anisotropy, and pion 
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intcrfcromctry in Au+Au collisions at ^/snn = 9.2 
GeV. The results are obtained from only about 3000 
events from the lowest beam energy run to date at 
the RHIC facility. The transverse momentum spec- 
tra of pious, kaons, and protons are presented for 
0-10%, 10-30%, 30-60%, and 0-60% coUision cen- 
trality classes. The bulk properties arc studied by 
measuring the identified hadron dN/dy, (pt), parti- 
cle ratios, vi (also at forward rapidity), W2, and HBT 
radii (i?out, Rside, and i?iong)- AH measurements are 
consistent with corresponding previous results from 
fixed target experiments at similar ^/snn- 

The (pt) for protons is higher than that for pi- 
ous, indicating some degree of collective flow in the 
radial direction. However, the difference between 
(Pt) for protons and kaons is considerably smaller 
at ^/SNN — 9.2 GcV than at ^/sNN = 62.4 and 200 
GeV at RHIC. This suggests that the average col- 
lective velocity in the radial direction at the lower 
beam energy is smaller compared to 62.4 and 200 
GeV collisions. 

The p/p ratio at midrapidity for ^sjvat = 9.2 
GcV collisions is much smaller, with a value of 
0.010 ± 0.001 (stat.) ± 0.003 (sys.), and the p/Tr+ 
ratio is larger compared to Au-|-Au collisions at 
y^SNN — 200 GcV. These measurements indicate 
large net-proton density at midrapidity in collisions 
at y/SNN = 9.2 GeV. In this region of high net- 
baryon density for 9.2 GeV collisions, the dominant 
channel for kaon production is the associated pro- 
duction. The /K^ ratio has a value of 0.38 ± 
0.05 (stat.) ± 0.09 (sys.) and the K+ /Tr+ ratio 
is slightly higher compared to that in collisions at 
^/i]^ = 200 GcV. 

The directed flow measurements, plotted as a 
function of pseudorapidity scaled by the beam ra- 
pidity, have similar values for three collision energies 
(Vs/vjv = 9.2, 62.4, and 200 GeV). A large vi sig- 
nal is observed at forward rapidities at y^s^N = 9.2 
GcV. These collisions could have significant contri- 
bution from protons that dominate at large \r]\ (spec- 
tator effects). The V2 measurements for charged 
hadrons, pious, and protons are also presented for 
= 9.2 GeV Au-hAu collisions at RHIC. The 
charged pion V2 as a function of pr is observed to 
be comparable with that from NA49 at similar col- 
lision energy. The STAR data at y^sjvTv = 9.2 GeV 
arc also found to follow the existing beam energy 
dependence of V2 for charged hadrons. 



The pion intcrfcromctry results give information 
of the size of the homogeneity region of the source. 
The pion HBT radii i?out, ^sido, and i?iong have val- 
ues 5.05 ± 0.96 fm, 3.52 ± 0.56 fm, and 3.25 ± 0.86 
fm, respectively. 

The kinetic freeze-out parameters are extracted 
from a blast- wave model fit to pion, kaon, and proton 
Pt spectra. We obtain Tkin = 105 ± 10 (stat.) ± 
16 (sys.) MeV and {(3t) = 0.46c ± 0.01c (stat.) ± 
0.04c (sys.). The chemical freeze-out parameters are 
extracted from a thermal model fit to the particle 
ratios at midrapidity. We extract T^h = 151 ± 2 
(stat.) ± 7 (sys.) MeV and pB = 354 ± 7 (stat.) ± 
30 (sys.) MeV for 0-10% central Au+Au coUisions 
at ./JWn = 9.2 GeV. 

These results from the lowest energy collisions 
studied up to now at RHIC demonstrate the capa- 
bilities of the STAR detector to pursue the proposed 
Beam Energy Scan. Large and uniform acceptance 
for all beam energies in a collider set up, excellent 
particle identification (augmented by the inclusion of 
a full barrel Time-Of-Flight [67] in addition to the 
large acceptance TPC), and higher statistics will of- 
fer significant quantitative and qualitative improve- 
ment over existing data. The QCD critical point 
program at RHIC will allow us to extensively ex- 
plore the QCD phase diagram. It will also allow us 
to search for the onset of various observations related 
to partonic matter that have already been uncovered 
at the highest RHIC energies. 
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